% 主程序：光学散射计算与分析
clear all;
close all;

% 定义基本参数结构体
params = struct();
params.theta1 = deg2rad(8);             % 入射角 (弧度)
params.lambda = 1319e-9;                % 激光波长 (米)
params.sigma2 = 10e-9;                  % 随机分量的粗糙度 (米)
params.A = 10e-9;                       % 周期分量的振幅 (米)
params.L = 10e-3;                       % 激光束尺寸 (米)
params.T = 10e-6;                       % 相关长度 (米)
params.n = 758;                         % 影响弓形的参数
params.m_max = 100;                     % 级数截断项数

% 验证参数条件
validate_parameters(params);

% 计算散射系数
[theta2_range, rho_rho_star, yy] = calculate_scattering(params);

% 分析结果
analyze_results(theta2_range, rho_rho_star);

% 绘制增强的结果图
plot_enhanced_results(theta2_range, yy);

% 进行参数敏感性分析
sensitivity_analysis(params); 